Adiabatic Mach-Zehnder interferometer via an array of trapped ions 



(N 

O 

(N 



(N 

a- 



(N 

>: 



X: 



Y. M. Feng\t and C. Lee^* 

^ State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, 
and Wuhan National Laboratory for Optoelectronics, Wuhan 430071, China 
^Graduate School of the Chinese Academy of Sciences, Beijing 100049, China and 
State Key Laboratory of Optoelectronic Materials and Technologies, 
School of Physics and Engineering, Sun Yat-Sen University, Guangzhou 510275, China 

(Dated: March 23, 2012) 

We explore the possibihty of implementing a Heisenberg-limited Mach-Zehnder interferometry via 
an array of trapped ions, which obey a quantum Ising model within a transverse field. Based upon 
adiabatic processes of increasing the Ising interaction and then decreasing the transverse field, we 
demonstrate a perfect transition from paramagnetism to ferromagnetic states, which can be used 
as the beam splitter for the multi-ion Mach-Zehnder interferometry. The achieved NOON state 
of the ions enables the Heisenberg-limited interferometry. Using currently available techniques for 
ultracold ions, we discuss the experimental feasibility of our scheme with global operations. 
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Since the foundation of quantum theory, as a natural 
result of superposition principle, quantum interference 
has attracted continuous attentions in both theoretical 
and experimental studies. It has been widely used to 
implement high-precision measurement, quantum infor- 
mation processing and so on. A well-known scheme for 
performing quantum interferometry is the Mach-Zehnder 
(MZ) interferometry, which has a beam splitter for split- 
ting the input states and another beam splitter for re- 
combining the output states. Up to now, quantum MZ 
interferometry has been accomplished via photons [1], 
electrons [2] , superconducting flux qubits [3] and trapped 
ions [4, 5] etc. 

Beyond the conventional quantum interferometry via 
unentangled states, it has been demonstrated that 
the measurement precision can be enhanced from the 
standard quantum limit (or the shot noise limit) to 
the Heiscnberg limit by using multipartite entangled 
states [6]. An excellent candidate is the NOON state 

i\N)a |0)b + |0)a 1^)6) /v^' ■^hich is an equal-probability 
superposition of all N particles in one of two paths 
denoted by a and b. The entangled ions for high- 
precision metrology have been proposed theoretically [7] 
and demonstrated experimentally [8, 9]. However, these 
schemes are subject to limited numbers of ions or the re- 
quirement for individual addressing. For an ensemble of 
thousands of neutral atoms, the possibility of performing 
a Heisenberg-limited MZ interferometry has been demon- 
strated via a quantized Bose-Josephson junction [10]. 

In this article, based upon the adiabatic processes and 
global operations on an array of ultracold ions, we present 
a realizable scheme for performing a Heisenberg-limited 
MZ interferometry. The ion array is described by a quan- 



tum Ising model of ferromagnetic (FM) interaction J 
and transverse field B, which has been used to simulate 
quantum magnetism [11-13] and quantum phase transi- 
tion (QPT) [14-16] by virtue of spin-dependent optical 
dipole forces [17]. In our scheme, if the system starts from 
a paramagnetic state dominated by B, we adiabatically 
increase J and then decrease B, and vice versa if the sys- 
tem starts from a FM state dominated by J. In contrast 
to tuning either J or i? in previous schemes [11, 16], the 
adiabatic processes in our scheme perfectly connects the 
two limits solely controlled by one of B and J. There- 
fore, theoretically, a pure NOON state can be prepared 
adiabatically from a SU(2) coherent state, which is the 
ground state for the system completely dominated by B. 
Our scheme requests only global operations, which is less 
challenging experimentally than individual addressing of 
the ions. By using the adiabatic process between param- 
agnetic and FM states as beam splitters, we are able to 
accomplish a Heisenberg-limited MZ interferometry via 
the NOON state of an array of trapped ions. 

Wc consider N ultracold ions confined in a linear Paul 
trap. Because only two ionic hyperfine states are used 
for implementing the MZ interferometry, each ion can be 
regarded as a spin-1/2 particle of two spin states \\.) and 
]t)- By globally addressing all the ions with particular 
lasers, the system obeys a transverse-field quantum Ising 
Hamiltonian [14-17], 



H 



i<j 



(1) 
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where cr* ^ Pauli operators for the i-th ion, B is the 

transverse field, and Jy — J/ \i ~ jf^ is the effective Ising 
interaction between ions i and j with J > denoting 
the nearest-neighboring interaction. For convenience, we 
consider the dimensionless model in units of S = 1 before 
the discussions of experimental possibility. This model 
has been experimentally realized by using either longi- 
tudinal [11, 17] or transverse modes [12-16] of the ions. 
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Obviously, there are two extreme cases: J = and B = 0. 
If J = 0, the ground state is a paramagnetic state of aU 
spins ahgned with the magnetic field, i.e., |-^^ • • • 
for B > or • • • ^) for B < 0. This ground state 

is a spin coherent state. If B = 0, there are two degen- 
erate ground states of all spins in either the spin-down 
state lii • • • i) or the spin- up state |tt • • ' t)- The equal- 
probability superposition of these two degenerate ground 
states is a NOON state. 
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exactly the FM states and the initial state is not the 
ground state for the second scheme [16]. 

To improve the fidelity to the desired NOON state, 
different from the single-step schemes of sweeping B or 
J, we adopt in the present work a two-step scheme ex- 
actly connecting two extreme cases, i.e., ( J = 0, B > 0) 
and (J > 0, B = 0), in which J and B are alternately 
changed in one of the two steps. Thus, the ground state 
for the initial Hamiltonian Hi = ^B'^^al. is exactly 
the spin coherent state along x— axis, and the NOON 
state is exactly a ground state for the final Hamilto- 
nian Hf = ^J2i<j Jij'^z'^i- the first step, the ratio 
J/ B increases from to 1 with B remaining unchanged. 
While in the second step, the ratio B/J decreases from 
1 to with J remaining unchanged. In our calcula- 
tions, the initial state is chosen as |— >— > ■ ■ ■ — >■), which 
can be easily prepared. The system first evolves with 
B = 1 and J ~ t/r from t = to t = t and then 
evolves with J ~ 1 and B = 2 — t/r from i = r to 
t — 2t. Therefore, for a sufficiently large r, the system 
will adiabatically evolve into the desired NOON state 
I*) = •••;> + Itt ■■ • t))/^/2 with very high fideli- 
ties, see Fig. 1(c) and Fig. 2. Besides, Fig. 1(d) shows 
that the NOON state achieved by our two-step scheme 
is much better than the ones achieved by the single-step 
schemes. 



FIG. 1: (Color online) Population evolutions of FM states 
for systems of different A^, in which the time t £ [0, 2t] with 
r = 5. (a) J = t/r and B = 0.4. (b) J = 0.4 and B = 2-f/r. 
(c) J = t/T and B = Hit € [0,r], and J = 1 and B = 
if t G [r, 2r] . (d) Population evolutions for N = 8 under 
different sweeping schemes. In (a), (b) and (c), the insets 
stand for the temporal modulations of B and/or J. 

To prepare the NOON state for the multi-ion MZ inter- 
ferometry under global operations, one may use a quan- 
tum adiabatic evolution connecting the initial spin co- 
herent state and the NOON state. Obviously, the initial 
state is the ground state for the system of ,/ = and 
B > and the desired NOON state is the ground state for 
the system of J > and B = 0. Two sweeping schemes 
have been suggested for preparing the NOON state with 
high fidelities. In the first scheme, J is adiabatically 
increased from to i?/J <C 1 (with constant B) [11]. 
In the second scheme, B is adiabatically decreased from 
B/J^ 1 to (with constant J) [16]. Starting from 
]—)■—> • • • — >), we show the dynamical populations in the 
two FM states \ii ■ ■ ■ i) and [tt ' • ' t), see Fig. 1(a) and 
(b). The two FM states always have the same popula- 
tions due to the absence of longitudinal fields and the 
fidelity to the desired NOON state is associated with the 
sum of these two populations. The final populations de- 
pend on the sweeping rate and the system size. But even 
if the time-evolution is perfectly adiabatic, for a finite- 
size system, the final fidelity to the NOON state can only 
approach unity but not exactly unity. This is because 
that the final ground state of the first scheme [11] is not 
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FIG. 2: (Color online) Fidelity to the NOON state versus the 
total sweeping time 2t for systems of different TV. Exper- 
imentally the fidelity is relevant to the ferromagnetic order 
parameter, which could be obtained by measuring the popu- 
lations of the spin-up and -down states with spin-dependent 
fluorescence [16]. 

On the other hand, at the expense of a longer total 
sweeping time (2t), the single-step scheme of tuning B or 
J may also achieve a desired fidelity. However, this is not 
experimentally practical for multi-ion systems restricted 
by the limited coherent time. Therefore, our two-step 
scheme of a shorter total sweeping time seems a better 
choice for preparing the desired NOON state. 

Utilizing the prepared NOON state, we may achieve a 
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Hcisenberg- limited MZ intcrfcromctry by using |tt ' ■ ■ t) 
and lii • • • i) as two paths and by global operations for 
all the ions. The schematic diagram for our adiabatic MZ 
intcrferometry is shown in Fig. 3, in which a free evolu- 
tion is sandwiched by two beam splitters achieved by adi- 
abatic processes connecting paramagnetic and FM states. 
The first beam splitter (BSl) is the adiabatic prepa- 
ration of the NOON state 1 5*) discussed above. Then 



the state 1^1') evolves into 



-iN4>/2 



in • • • t> + 



e'^-^/^ III . . . 1))/^ in the following free evolution of 
time duration T under the government of Hq = cjq '^1 ^ 
where wq is the transition frequency between \'f) and |^) 
and the accumulated phase </> = ujoT- Lastly, the sec- 
ond beam splitter (BS2) is accomplished by the inverse 
process of the BSl. 



■>> BSl 



,\n.A) 

\u...i) 



Detector 



=0 



FIG. 3: (Color online) The schematic diagram for the adia- 
batic MZ interferometry via trapped ions. Here, BSl prepares 
the NOON state via adiabatic processes connecting paramag- 
netic and FM states, then two paths accumulate a relative 
phase in a following free evolution, and BS2 recombines the 
two paths for interference. 



To extract the accumulated relative phase between the 
two paths, one has to transform the phase information 
into the amplitude information of the final state itself 
or the expectation information of a particular observable 
for the final state. For systeins of few trapped ions, this 
could be done by a controUed-NOT gate with the first 
ion being the control qubit and the rest being the tar- 
get ones followed by a Hadamard operation on the first 
ion [7]. However, the individual addressing used in this 
procedure is very difficult to be accomplished in systems 
of large numbers of particles. In our scheme with full 
global operations, we may accomplish this procedure in 
another adiabatic process, which is the inverse process of 
the BSl, i.e., by first slowly increasing B/J front to 1 
(with J fixed) and then slowly decreasing J/ B from 1 to 
(with B fixed). As a result, for extracting the phase 
shift from our multi-ion MZ interferometry, we have to 
measure populations of the ground state |0) and the first 
excited state |1) of the Hamiltonian Hi [18]. Under a 
global Hadamard operation, |0) turns to be spin- up with 
all the ions, and |1) becomes an entangled state with 
one spin-down and others spin-up. By applying a probe 
laser for coupling | \.) to an excited level, it is possible 
to distinguish |1) from |0) by detecting spin-dependent 
fluorescence signals [16]. Therefore, we could obtain the 
population of the first excited state. 



and the corresponding phase sensitivity is 

= AP,/{d(p,) /a0) = 1/N, (3) 



Pi = sin2(iV0/2), 



(2) 



with the population operator Pi = |1) (1| and its variance 
APi. This is the Hcisenberg limit, where the relative 
phase (j) is measured inore precisely than the disentangle- 
inent case by \/N times. Since (p = (jJqT and ojq depends 
on the atomic configuration and the local magnetic field, 
our Hcisenberg-limited MZ interferometry is of practical 
applications in understanding the atomic configuration 
and measuring the local magnetic field. 

Experimentally, the transverse-field quantum Ising 
inodcl could be realized by an array of ultracold ions 
which are globally illuininated with off-resonant lasers 
and resonant Raman beams [12-17]. Because of the off- 
resonant lasers, the inter-ion Ising interactions are in- 
duced by the spin-dependent forces [19] with the assis- 
tance of phonon modes. The effective transverse fields 
arc generated by the resonant Raman beams, which res- 
onantly couple the two spin states. In the rotating frame, 
the ionic array under the Lainb-Dicke limit is equivalent 
to a quantuin Ising system within a transverse magnetic 
field. With the control of the frequencies and intensi- 
ties of the off-resonant lasers, the sign and strength of 
the Ising interaction can be dynamically adjusted [15]. 
By tuning the strength of the resonant Raman beams, 
we can dynamically adjust the strength of the transverse 
field [13, 14, 16]. 

For an array of ^^^Yb+ ions, two clock states 
\F = 0,mF = 0) and \F = 1, = 0) are em- 
ployed as the two spin states \l) and |t) for a spin-1/2 
particle, respectively. Based upon an array of ultracold 
171 Yb -I- ions in a linear trap, the QPT in the transverse- 
field quantum Ising model has been simulated [16], in 
which both B and J may be larger than the order of 
kHz. Therefore, we assume the maximum value for B 
and J is Bq = Jq = 50 kHz for our two-step sweep- 
ing scheine shown in the inset of Fig. 1(c). To coinpare 
with the dimensionless model, we transform the original 
Hamiltonian H into H/ Jq, which means that the energy 
is in units of Jq and the time is in units of I/Jq. In the 
first step, wc fix _B = 50 kHz and slowly vary J according 
to J = Joxt/r from t ~ to t = t, where r — 100/^s. 
Then we fix J = 50 kHz and slowly sweep B following 
B = BqX (2— t/r) from t~Ttot = 2r. Besides, to main- 
tain the coherence of the NOON state, the free evolution 
time T should be much shorter than r. Therefore the re- 
quired operations from BSl to BS2 in our proposed MZ 
interferometry take about 400 /iS, which is shorter than 
the gating time (longer than 500 /is) in Rcf. [16] and 
should be feasible with currently available techniques for 
finite spin size. 

Our scheme could also be realized by using other types 
of ions. For an array of ^'^Ca"'' ions in a linear trap, | l) 
and I t) are encoded by two hyperfine states 1 5*1/2, ~ 
-1/2) and \Si/2,ms = 1/2) [or [5i/2, mg = -1/2) and 
the mctastable state \Dc,/2, ms ~ —1/2)], respectively. 
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It has demonstrated the entanglement of fourteen Ca"*" 
ions in a linear trap [20]. More recently, the universal 
digital quantum simulation based upon a transverse-field 
quantum Ising model has been implemented by using 
^"Ca^ [21]. For a system of fourteen '*°Ca~'" ions, our 
intcrfcrometry scheme could enhance the measurement 
precision by nearly four times in comparison to the stan- 
dard quantum limit. 



(a) (b) 
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FIG. 4: (Color online) The undesirable impact from the de- 
viation in the process of BSl for eight ions by sweeping time 
t = 2r = 200^s. (a) Populations in |tt ' ' • t) and |i| ■ • ■ i) 
with respect to 5/ Jo, and (b) the fidelity of the prepared state 
to the ideal NOON state with respect to S/Jq with the vari- 
ation of B and J in the inset during the adiabatic evolution. 

Besides unpredictable imperfection in operations, the 
main errors in realistic experiments are caused by spon- 
taneous emission and intensity/frequency fluctuations of 
the Raman beams, which lead to decoherence and sup- 
press the populations of FM states [16]. The spontaneous 
emission in our scheme can be minimized by enlarging 
the Raman detuning from the excited state, provided 
that we increase the laser power accordingly to main- 
tain the original Ising interaction and the transverse field 
strength. The intensity fluctuations of the Raman beams 
only yield AC Stark shift with 2% rms error [12], which 
is not serious in our scheme. Alternatively, we may em- 
ploy *°Ca+ ions irradiated by a 729 nm laser to couple 
the ground state lS'i/2i ™s = ~l/2) with the metastablc 
state lI?5/2, ms ~ —1/2), which completely remove the 
necessity of Raman beams in manipulation [22]. For 



the operational imperfection, we simply consider a typ- 
ical one in the non-resonant carrier transition when we 
tunc the parameter B, which yields unwanted longitu- 
dinal bias field S al and may cause a deviation from 
the expected NOON state. This imperfection is shown 
in Fig. 4, and it could be suppressed by carefully adjust- 
ing the direction and the polarization of the laser beams 
and/or by using additional pulses [11] to compensate the 
deviation. 

In our scheme, the two beam splitters form a loop- 
like adiabatic operation, which may generate a geometric 
phase [23] . Fortunately, the generated geometric phase is 
a constant relevant to the loop area, but do not depend 
on the free evolution time T (i.e. the accumulated phase 
luqT) [24]. Therefore, by measuring the relative phase for 
different free evolution time, we may eliminate the infiu- 
ence of the geometric phase by reducing a common shift 
from measurement results and the measurement precision 
of ujq is thereby independent from the geometric phase 
accumulated in the adiabatic operation. 

In conclusion, we have proposed a simple and practi- 
cal scheme to carry out a Heisenberg-limited MZ inter- 
ferometer with an array of trapped ions, which obey a 
transverse-field quantum Ising model. The multi-ion MZ 
interferometry is implemented by only global operations, 
which favors scalability. However, more challenges would 
appear if more ions are involved, such as less homoge- 
neous laser irradiation on the ions, smaller energy gap in 
the Ising model and weaker laser-ion coupling. Neverthe- 
less, for finite ions, the global operation is less difficult 
experimentally than individual addressing. Provided the 
nearly perfect adiabaticity, our scheme includes a QPT in 
each beam splitter. By suppressing unpredictable errors 
in realistic experiments and elaborately modifying the 
results, our scheme for quantum metrology is in princi- 
ple able to reach the ultimate precision limit beyond the 
standard quantum limit [25]. 

This work is supported by the NBRPC under Grants 
No. 2012CB821305 and No. 2012CB922102, the NNSFC 
under Grants No. 10974225, No. 11075223 and No. 
11004226, the NCETPC under Grant No. NCET-10- 
0850 and the Fundamental Research Funds for Central 
Universities of China. 



[1] J. G. Rarity, P. R. Tapster, E. Jakeman, T. Larchuk, R. 

A. Campos, M. C. Teich, and B. E. A. Saleh, Phys. Rev. 

Lett. 65, 1348 (1990). 
[2] Y. Ji, Y. Chung, D. Sprinzak, M. Heiblum, D. Mahalu, 

and H. Shtrikman, Nature (London) 422, 415 (2003). 
[3] W. D. Oliver, Y. Yu, J. C. Lee, K. K. Berggren, L. S. 

Levitov, and T. P. Orlando, Science 310, 1653 (2005). 
[4] R. Huesmann, Ch. Balzer, Ph. Courteille, W. Neuhauser, 

and P. E. Toschek, Phys. Rev. Lett. 82, 1611 (1999). 
[5] D. Leibfried, B. DeMarco, V. Meyer, M. Rowe, A. Ben- 

Kish, J. Britton, W. M. Itano, B. Jelenkovic C. Langer, 

T. Rosenband, and D. J. Wineland, Phys. Rev. Lett. 89, 



247901 (2002). 

[6] V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 

1330 (2004); Phys. Rev. Lett. 96, 010401 (2006). 
[7] J. J. Bollinger, W. M. Itano, D. J. Wineland, and D. J. 

Heinzen, Phys. Rev. A 54, R4649 (1996); S. Boixo, A. 

Datta, M. J. Davis, A. Shaji, A. B. Tacla, and C. M. 

Caves, Phys. Rev. A 80, 032103 (2009). 
[8] C. A. Sackett, D. Kielpinski, B. E. King, C. Langer, V. 

Meyer, C. J. Myatt, M. Rowe, Q. A. Turchette, W. M. 

Itano, D. J. Wineland, and C. Monroe, Nature (London) 

404, 256 (2000). 
[9] D. Leibfried, M. D. Barrett, T. Schaetz, J. Britton, J. 



5 



Chiaverini, W. M. Itano, J. D. Jost, C. Langer, and D. 
J. Wineland, Science 304, 1476 (2004). 
[10] C. Lee, Phys. Rev. Lett. 97, 150402 (2006); C. Lee, J. 
Huang, H. Deng, H. Dai, J. Xu, Front. Phys. 7, 109 
(2012). 

[11] A. Friedenauer, H. Schmitz, J. T. Glueckert, D. Porras, 
and T. Schaetz, Nature Phys. 4, 757 (2008). 

[12] K. Kim, M.-S. Chang, R. Islam, S. Korenblit, L.-M. 
Duan, and C. Monroe, Phys. Rev. Lett. 103, 120502 
(2009). 

[13] K. Kim, M.-S. Chang, S. Korenbht, R. Islam, E. E. Ed- 
wards, J. K. Freericks, G.-D. Lin, L.-M. Duan, and C. 
Monroe, Nature (London) 465, 590 (2010). 

[14] E. E. Edwards, S. Korenbht, K. Kim, R. Islam, M.-S. 
Chang, J. K. Freericks, G.-D. Lin, L.-M. Duan, and C. 
Monroe, Phys. Rev. B 82, 060412(R) (2010). 

[15] G.-D. Lin, C. Monroe, and L.-M. Duan, Phys. Rev. Lett. 
106, 230402 (2011). 

[16] R. Islam, E.E. Edwards, K. Kim, S. Korenblit, C. Noh, H. 
Carmichael, G.-D. Lin, L.-M. Duan, C.-C. Joseph Wang, 
J. K. Freericks, and C. Monroe, Nat. Commun. 2, 377 
(2011). 

[17] D. Porras and J. I. Cirac, Phys. Rev. Lett. 92, 207901 
(2004); X.-L. Deng, D. Porras, and J. I. Cirac, Phys. Rev. 
A 72, 063407 (2005). 

[18] For N ions, the ground state |0) = [([ t) + i ;))/%/2]®^ 
and the first excited state |1} consists of N degenerate 
sub-states. For a system of N = 3, its three degenerate 
sub-states are |ei) = (-] m) + \ til)-! ItD + l m))/2, 



|e2) = (-1 tti) - I m) + I m) + I iit))/2 and jea) = 
(-1 t;t) - i m) + I m) + \ W»/2. under a global 
Hadamard gate, |0} turns to be j t)®^ and each sub-state 
of |1) is an entanglement of two spin-up and one spin- 
down. We may accomplish the readout by detecting the 
spin-dependent spontaneously emitting photons. 

[19] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. 
Mod. Phys. 75, 281 (2003). 

[20] T. Monz, P. Schindler, J. T. Barreiro, M. Chwalla, D. 
Nigg, W. A. Coish, M. Harlander, W. Hansel, M. Hen- 
nrich, and R. Blatt, Phys. Rev. Lett. 106, 130506 (2011). 

[21] B. P. Lanyon, C. Hempel, D. Nigg, M. Miiller, R. Ger- 
ritsma, F. Zahringer, P. Schindler, J. T. Barreiro, M. 
Rambach, G. Kirchmair, M. Hennrich, P. ZoUer, R. Blatt, 
and C. F. Roos, Science 334, 57 (2011); D. P. DiVin- 
cenzo. Science 334, 50 (2011). 

[22] M. Riebe, K. Kim, P. Schindler, T. Monz, P. O. Schmidt, 
T. K. Korber, W. Hansel, H. Haffner, C. F. Roos, and R. 
Blatt, Phys. Rev. Lett. 97, 220407 (2006). 

[23] M. V. Berry, Proc. Roy. Soc. Lond, Ser. A 392, 45 (1984); 
F. Wilczek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984); 
Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 
(1987). 

[24] E. Sjoqvist, A. K. Rati, A. Ekert, J. S. Anandan, M. 

Ericsson, D. K. L. Oi, and V. Vedral, Phys. Rev. Lett. 

85, 2845 (2000). 
[25] D. J. Wineland and D. Leibfried, Laser Phys. Lett. 8, 

175 (2011). 



